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Abstract 

A numerical approach to compute tensor integrals in one-loop calculations is presented. 
The algorithm is based on a recursion relation which allows to express high rank tensor 
integrals as a function of lower rank ones. At each level of iteration only inverse square 
roots of Gram determinants appear. For the phase-space regions where Gram determi- 
nants are so small that numerical problems are expected, we give general prescriptions 
on how to construct reliable approximations to the exact result without performing 
Taylor expansions. Working in 4 + e dimensions does not require an analytic separa- 
tion of ultraviolet and infrared/collinear divergences, and, apart from trivial integrals 
that we compute explicitly, no additional ones besides the standard set of scalar one- 
loop integrals are needed. 
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1 Introduction 



Computing radiative corrections in high-energy Particle Physics demands an increasing ca- 
pability of manipulating complicated objects. When the number of particles undergoing the 
scattering process grows, formidable complications arise. For example, to date no complete 
QCD or Electroweak (EW) calculation exists at the one-loop level describing processes in- 
volving 2 ingoing and 4 outgoing particles. On the other hand, multi-particle processes need 
to be studied with great accuracy at the next generation of pp and e + e~ colliders, hence 
we cannot escape the technical subject of efficiently compute radiative corrections, at least 
at the one-loop level. To do this three main problems must be faced: the large number of 
Feynman diagrams, the reduction of tensor integrals to scalar ones, and the control over the 
numerical inaccuracies. 

As for the first issue, while several tree level algorithms to compute amplitudes without 
making explicit reference to Feynman diagrams exist pQ , so far there is no equivalent working 
technique when loops enter the game. At any rate, for moderate values of external particles 
the number of contributing Feynman diagrams can be in principle manageable. For example, 
this number is of the order of thousand for the EW process e + e~ — > fi'u^ud, depending on 
the chosen gauge. What really renders the calculation difficult is the fact that each diagram 
still requires a lot of work to be computed. This is the second issue mentioned above. 

In order to be concrete, let us sketch out a typical one-loop calculation performed in 
n = 4 + e dimensions. The corresponding amplitude can be written as a sum of tensor 
integrals 7^ times external tensors S. 

A(n) _ ft™) CW'W fl ^ 

m,j,i 

where 

4S;,x- w = f dn <l nn ' 9 n > D k = (q + p k f - k = 0,--,m, p£ = 0, (2) 

with j labelling the different momenta p k entering and masses m k running the loop, and where 
St^J^ only depend on the external kinematics. The (m + l)-point tensor one-loop integrals 
are usually evaluated in two steps. First one writes down the most general decomposition 
in terms of the external momenta and the metric tensor 

c; :/ „.„ e <f (s,s ir ..,s m ){[gr\ Pl r...\p m r-}^ lH , (3) 

2s + J2s k = i 

where {[g] s [pi] si • • • [Pm] Sm } Atl ... At . corresponds to the symmetrical tensor combination, each 
term of which is constructed from s metric tensors, Si momenta pi, ■ ■ •, and s m momenta 
p m . For example, 

= flVnwPi^ + 9mnzPx m + ■ ( 4 ) 
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Secondly one looks for explicit expressions for the scalar coefficients cf (s, S\, ■ ■ ■ , s m ). In 
the Passarino-Veltman treatment j2] they are expressed in terms of a minimal set of scalar 
one-loop integrals with denominators raised to the first power only [3] 

In other methods |U E] the decomposition is performed in terms of an enlarged set of scalar 
integrals with shifted dimensionality and denominators raised to generic powers 

j{n+v) = f An+u) 1 /gX 

Of course, there are relations among the integrals in Eq. (|BJ) and the minimal set in Eq. 
(0), which have been extensively studied in Ref. [I]. For large m values the algebraical 
complexity of the described methods quickly becomes overwhelming, mainly due to the 
rapidly increasing number of kinematic variables present in the problem. For example, the 
application of Eq. (jHJ) to 1^1 ^u,^ generates 24 scalar coefficients. 

An alternative to this procedure is a numerical approach. The ideal situation would be 
that one simply writes down Eq. while all the rest is handled by a numerical program. 
The main problem with such an approach is handling ultraviolet and infrared/collinear 
singularities. Along this road a complete formalism has been recently presented in the 
framework of QCD by Giele and Glover In their method one first analytically separates 
the divergent contributions arising from the tensor integrals in Eq. (j2J), using the techniques 
of Refs. [Z1IB1IH]- Then one computes numerically the kinematical coefficients of the resulting 
finite 4-dimensional integrals. A different approach is sewing tree amplitudes together to 
construct loop amplitudes, as proposed in Ref. [TO]. Another way is combining virtual 
and real contributions to cancel the divergences in the loop integration [TI], or constructing 
counter-terms diagram by diagram ^2]- Finally, the authors of Ref. ^3] develop a pure 
numerical approach where in the Feynman parameter space any one-loop integrand is cast 
in a form well suited for numerical calculation, in such a way that all possible divergences 
get automatically extracted. 

In this paper we present a new method in which almost all the work can be performed 
numerically: the tensor integrals in Eq. (J2J) are numerically reduced to the minimal set in Eq. 
(JHJ), and the ultraviolet and infrared/collinear divergences are controlled without performing 
any explicit subtraction. Then, any amplitude can be calculated simply contracting the 
numerically computed tensor integrals with the external tensors S m ,j (see Eq. Our 
main result will be the derivation of the set of recursion relations that link high rank tensor 
integrals to lower rank ones, and which allow the numerical reduction of the former to 
the minimal set in Eq. (j3J). On the other hand, our solution to the problem of handling 
divergences is splitting beforehand any tensor into its pure 4-dimensional part plus any other 
additional contributions, which are trivial to compute. Such a procedure renders unnecessary 
an analytical separation of the divergent parts, as required by the recursion relations of Ref. 
[7j, and minimizes the analytic work. At the end all divergences are contained in the pole 
parts of the scalar integrals, parts to which one can give any value to numerically check 
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all relevant cancellations. Furthermore, internal and/or external masses do not pose any 
particular problem, so that the method can be applied to both QCD and EW calculations. 

The tensor integrals are usually well behaved when two or more momenta become lin- 
early dependent, as observed in Refs. [T H I15 | ITT)]. However, in general reduction formulae 
introduce inverse powers of Gram determinants 



in the decomposition terms. So when they go to zero, large cancellations among the different 
terms must take place, giving rise to numerical instabilities. This is the third problem 
one has to face. In Ref. J7| this is solved by systematically building up combinations of 
well behaved functions in the limit of vanishing Gram determinants. The drawback of this 
approach is that the class of needed loop functions has to be enlarged, including integrals 
computed in higher dimensions, such as those in Eq. (jHJ). When things are re-expressed 
in terms of (4+e)-dimensional integrals, Gram determinants are reintroduced and explicit 
Taylor expansions are required to deal with the problematic phase-space regions. A possible 
solution is presented in Ref. |18j . where an extrapolation from the inner phase-space region 
is used for the dangerous points. In the numerical approach of Ref. ^3] all scalar coefficients 
c w j n Eq. (j2J) are cast in a form well suited for numerical evaluation. In Ref. ^H] 4- 
dimensional pentagon like tensors are reduced to box like tensors, avoiding the occurrence 
of rank m = 4 Gram determinants. Finally, the case of exactly zero Gram determinant is 
solved in Ref. [7j using the pseudo- inverse of the Gram matrix. 

In our method part of the Gram determinant singularities compensates in such a way 
that at each step of the iteration only inverse square roots of Gram determinants appear, 
improving the numerical stability of the calculation. Moreover, the expressions can be very 
naturally arranged in groups which are well behaved in the limit of linearly dependent ex- 
ternal momenta, allowing to keep as local as possible the numerical cancellations that occur 
among the loop functions. The case of exactly zero Gram determinants can be treated with- 
out any problem, and for the problematic regions of nearly vanishing Gram determinants, 
we show how to systematically construct reliable approximations to the exact result. Unlike 
in Ref. building up such approximations does not require explicit Taylor expansions. 

The paper is organised as follows. Our master formula for the general 4-dimensional case 
is presented in next Section. The algorithm is extended to the n-dimensional case in Section 
El The particular cases of 3-point and rank one tensors, to which our general formula does 
not apply, are discussed in Sections 0] and respectively. Whereas in Section El we perform 
a detailed study of the dangerous collinear and coplanar configurations. Section [7| is devoted 
to conclusions; and technical details are worked out in three Appendices. 

2 The 4-dimensional method 

In the following we shall show our master formula to reduce high rank tensor integrals to 
lower rank ones. Let us first clarify the notation. Throughout the paper we drop the index 
(n) from loop integrals when the space-time is 4-dimensional. We also omit the index j for 




1 < i, j ' < m 



(7) 
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simplicity, although it is understood that an (m + l)-point integral depends on m external 
momenta and m + 1 internal masses, m-point tensor integrals obtained from an (m + 1)- 
point one by eliminating the k th denominator are written specifying the index of the dropped 
denominator as an argument 

..... (*) = / *q— — ; (8) 

and as usual tensor indices are raised and lowered with the metric tensor. 

Our main result is the recursion relation for reducing 4-dimensional tensor loop integrals 
with m > 2 and rank higher than 1, 



Im; fivp—r — ^ TpvXct \^rn;p—r\ ^ Tp,v \wi§Im; p—r ~t~ Im—1; p"-r(0) } 

— J^^pi^A ^/ 30 ^m;p— t ^m-l;p— t(^) — ^m-l;p— t(^) ~£*3a <^ m; p... r | > 



where 



Jm;p-T — (fw r 2 + f20 r l)Im ]P - T + r 2^m-l;p-r(l) + r l p-r (2) 

- (^ + r 2 A )C_ 1;p ... T (0), (10) 

and T^xa, T^u\, T^v and r^ 2 only depend on three linearly independent external momenta, 
which we assume to be Pi,2,3 • The scalar factors (3 and 7 are only functions of pi )2 ; whereas 

/fco = m l~ m l-vl- (11) 

Eq. (jUJ) can then be iterated to compute numerically all 4-dimensional tensors I m ,j;ni-in> 
without explicit tensor decomposition, starting from the standard set of scalar loop functions 
in Eq. (JSJ). Notice that the variable shift q — > q — p± is needed before applying the next 
iteration to any tensor integral having (0) as an argument. 

In order to derive Eq. ©, we need to express the product of two loop momenta as a sum 
of terms with at most one loop momentum times loop denominators, or internal masses, and 
external tensors (Eq. (|24|) below). First, we write the internal momentum as a sum of four 
judiciously chosen massless vectors 

q = J2c i £ l , 4 2 = 0- (12) 
i=i 

Following Ref. [20] we construct £i j2 from two independent external momenta, which we 
assume to be pi >2 , 

Pi = £i + «i£ 2 , p 2 = £ 2 + a 2 ^i- (13) 
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Then 2 

4 = Pipi ~ OtlPl), 4 = P(P2 ~ "2 Pi), 

ai = ((Pi ■ p 2 ) T \/A)/pl A = -A 12 = (p! • p 2 ) 2 - p 2 p 2 , (14) 

«2 = a\ Pl/pl, 13 = 1/(1- a 1 a 2 ). 

Note that £1.2 can be also complex. The other two independent massless vectors 4,4 are 
taken to be 

^ = ^ 1 )fa)-u(4), ^ = C(^)fw-it(fi) l w- = (l-7 5 )/2, (15) 
thus fulfilling the equalities (see Appendix 1X1) 

(4,4 ■ 4,2) = *1 = tl = , (4 ■ 4) = -4(4 ■ £2) ■ (16) 

Using them the coefficients Cj in Eq. (J12)) can be simply written 

(g • 4) _ (g-_4) (g • 4) (g ■ 4) n7 , 

1 (4-4)' 2 (4-4)' 3 " 4(4-4)' 4 ~ 4(4-4)' 1 } 

Now, it is convenient to distinguish between £ 12 and 4,4 i n Eq. (|12|) for the contribution 
of the first two 4-vectors can be expressed as a sum of products of denominators, internal 
masses and external momenta, 

P n _J_ n 

7 2 1 

= i [2 (g •4)4 M + 2(g-4)4 M ] 

= {/10 r 2/, + /20 r i M + A r 2 M + A r lM - D (r^ + r 2At )} , 

= (?-4)4 / , + (g-4)4 M , 7 = 2 = 2(4.4), (18) 

1 + ai a 2 



where we have made use of 



(q-Pk) = ~[D k -D + f kQ \ , (19) 



When one or both p\ vanishes, ,9=1, and 



p 2 = 0, p|^0 =► ai = 0, «2 = 1^/(2 (pi- pa)), 

p?^0, p 2 2 = a 1 =p 2 1 /(2(p 1 -p 2 )) ) a 2 =0, 

pf = 0, p\ = => ai = 0, a 2 =0. 



These limits are smoothly approached taking for a\ the solution with -yA (+vA) when (pi ■ p 2 ) > 
((pi ' P2) < 0), which is what we stand for =py/~A. 
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with fko defined in Eq. (JUJ), and 

r 1 = (£ 1 -ai£ 3 ), r 2 = [£ 2 ~ a 2 £1) ■ 
Then applying Eq. ()18j) twice and symmetrizing on \w, we find 



(20) 



t 



L 
27 



D x q a 



20 



27 



(21) 



The factor in square brackets contains reconstructed denominators while the factor in curly 
brackets only depends on the external kinematics. The QpQ v term can be also decomposed 
using the properties of the 4- vectors £3,4 (see Appendix|X]for details) 

(q-h)(q-h) = 4(q-£ 1 )(q-£ 2 )-2q 2 (£ 1 -£ 2 ) = f3q a D a - 1 q 2 ) 

2 



(?•%)) (?•%)) 



(b-£ 



•*4(3). 



{[q 2 (£i-£2)-2{q-£ 1 )(q-£ 2 )}(b-£ m ) 



+ 2 [(g • 4)(& " £2) ~ (q ■ b)(h ■ £ 2 ) + (q ■ £ 2 )(h ■ b)] (q ■ 4(4))} 



[l(q-b)-(3b a D a ](q-£ m )} , V6^i, 2 . (22) 
Now, choosing the arbitrary 4- vector b = p s and reconstructing again denominators, we get 
Q,Q U = /3[DP(flig fn T ia ,}-'y[D + ml]{T ia ,} 



+ 



2(3p 3a [D a q x ]~ 1 (D 3 -D + f 30 )q x {T^ A } 



T, 



3/i<-3!^<-4A «4^<-4y<-3A 



(4-p 3 ) (4 -Pa) 



(23) 



Finally, inserting Eq. (|2*3j) into Eq. (|2*T|) we obtain the decomposition relation for the product 
of two loop momenta 



£ 
27 

1 

47 



D x q a 



T, 



1 



[luXa 



47 

(D 3 -D + / 30 ) 



A) + ™o 

2/5 
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P3«^ a )q X 



T, 



[J,uX ) 



where 



9fJ.X 9v 



27 



9Xa 



27 / 27 



/if 



(24) 



(25) 



Then Eq. © follows after dividing Eq. (|24|) by the (m + 1) denominators in Eq. (J2J), 
multiplying by the remaining 4- vectors q p - ■ ■ q T and integrating over d A q. As it stands, Eq. 
(JHJ) is valid in non exceptional phase-space regions where pi,2,3 are all independent, otherwise 
zeros occur in the denominators. In Sectional we shall consider the case of exceptional or 
nearly exceptional configurations. The reason why it is convenient to group terms in the 
way we have done will be clear there. 

As a last remark we observe that the rank of the tensor integral on the l.h.s. of Eq. 
should be at least 2 and m > 3, otherwise there are not enough external momenta to 
perform the decomposition neither denominators to cancel the reconstructed ones. These 
two particular cases are discussed in Sections [*"] and EJ respectively. 



3 The method in n dimensions 

The derivation of Eq. (""]} breaks down when working in n dimensions because Eq. f)12jl 
is only valid in 4 dimensions. To deal with the n-dimensional case we need to reduce the 
problem to 4-dimensional tensors, that we know how to handle. As in Refs. |2U[ |2*T| "2*2~| 12*3] 
only unobservable objects are considered to live in n = 4 + e dimensions, with 4 and e- 
dimensional quantities always orthogonal to each other. In particular, only the integration 
momentum q is n-dimensional in our integrals. For notational purposes from now on we put 
a bar over n-dimensional quantities and a tilde over e-dimensional objects. For example, 

q = q + q, q 2 = q 2 + q 2 , (26) 

where q is purely 4-dimensional. Being more explicit, the rank i n-dimensional tensor inte- 
grals we want to evaluate read 



= / d n q^ j- "Q % k , D k = (q + Vk f - m\ . (27) 

J Uq ■ ■ ■ JJ m k=1 

Using Eq. (j"""""j) to split the numerator momenta of Eq. (|27jl. we get 



n i n = n % k + e ?a^i n 

k=l k=l si=l k^si 

i i i 

+ e ^ %s 2 n % k + • • • + n %u ■ (28) 

s\>S2 fc^si,S2 k=l 

Since the momenta p k are purely 4-dimensional, all terms containing an odd number of q in 
the numerator vanish, and all terms with an even number of q can be only proportional to 
symmetric combinations of the metric tensor in e dimensions 

9/j.v = 9p,i> i 9fivp<j — 9p,v 9pa 4" 9p.p 9ua 4" 9fj.tr 9vp > etc. (29) 

Therefore 

= / ^ — 7T n ^. n ^ = . (3°) 
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and 

4^-™-^ = r (2) 2 'T{{le)I '■■ ■ !>,„ M ' (31) 

where q 2 = q^ and t = 0, 1, 2, ■ ■ -. As it is clear from the previous equation, a new class 
of integrals involving q 2i appears. We will use for these objects the notation j("- ;2£ ). For 
example, the last integral in Eq. (p?T]) is denoted by ^ ; ^f?.. M , • Such integrals are very easy 
to evaluate at Oil), and we do it in Appendix IB1 One can argue that, since a pole in 1/e 
appears in Eq. (J3"T]) . higher order terms might be needed. However, the e-dimensional tensor 
indices in Eq. (J31|) can only be contracted with e-dimensional indices (typically combinations 
of tensors), otherwise they give zero because of the orthogonality between 4-dimensional 
and e-dimensional spaces. In this contraction the pole cancels out, so that 0(e) terms can 
be safely neglected in the physical limit e — * 0. 

Finally, the purely 4-dimensional integral ij^^...^, coming from the first term of Eq. 
(|28|). can be reduced using Eq. @ with / — > I^ n > 3 . However, there is still one important 
modification we have to take care of. The 4-dimensional denominators appearing in Eq. (|24j) 
differ from the n-dimensional ones by an amount q 2 

D k = D k - q 2 . (32) 

To compensate for this, the only replacement needed in Eq. Q, before applying it to the 
n-dimensional case, is 

m 2 T ^rn 2 r(«) _ T(n;2) (oo\ 

l'l'0 1 m;p-T —> << l Q 1 m\p-T 1 m;p-r \°° ) 

Summarizing, any n-dimensional one-loop amplitude can be written 

A(n) _ \p An) _ (oa\ 
/ 1 m,j; pi---fii "raj ■ V / 

m,j,i 

After splitting the n-dimensional momenta according to Eq. ()28|). 1^1.^.. .„. can be evaluated 
with the help of Eq. (jHJ), together with the replacement given in Eq. (JHSj) . The additional 
integrals are always of the type 1^^...^- Since, after all, they are also 4-dimensional tensors, 
it would be in principle possible to compute them with the help of Eq. Q. In practice, the 
direct computation given in Appendix El is more convenient when one is interested in taking 
the limit e — ► 0. 

To conclude, we notice that the technique of splitting loop tensors in 4-dimensional plus 
e-dimensional parts can help also outside the algorithm we are presenting here. For example, 
the method to reduce pentagon tensor integrals to box tensor ones presented in Ref. JH] 
relies on 4-dimensional objects. However, a strategy for a possible extension to n-dimensions, 
which requires an explicit subtraction of soft and collinear divergences fHJ|21], is outlined 
by the authors. Instead, as described above, via the splitting in Eq. (|28|h the e-dimensional 

3 Of course, the tensors in the numerator of stay 4-dimensional. 
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part of the tensors can be computed separately; and the algorithm of Ref. [THj directly 
applied to the remaining 4-dimensional part of the tensor integrals (1^1^..^ in our notation) 
without any need of introducing a regulator A in the intermediate stages of the calculation. 
The only additional modification is the replacement in Eq. (|32jl for all the reconstructed 
denominators appearing in Ref. |19j . 



4 Three-point tensors 

Eq. (jHJ) cannot be applied when m = 2 because of the lack of a third 4-momentum p^ to 
reconstruct denominators. In this Section we derive a specific recursion relation for this case, 
which is valid for rank 2 and rank 3 three-point tensor integrals: 4 

An) _ ]3_ T , r T(n)A<n J_, r 2 An) , An) _ An; 2) "I /or^ 

J 2;Mp) - O^y [iv\o \^ 2; (p) J , V \ m J 2; (p) + J l; (p) l U J J 2; (p) / > 



;(p) 



where J^fl^ is the combination of one-loop integrals given in Eq. (fTUj) and 



T ^Act = ^flW - ^ j + 0„A (#p CT - ^ j + 7^ V • ( 36 ) 

Things get more complicated for higher rank tensors and the general solution is given in 
Appendix O 

To derive the former recursion relation, we shall make use of the following 
Theorem: 



D D 1 D 2 
1 

'd D 1 D 2 



d n q - i (q- £ 4 y = , Vi= 1,2,3- -•. (37) 



Proof: 



^D^D 2 {q ' = ® ' ' ' • (38) 

On the other hand, the tensor integral I^ui-m acUT iits a decomposition in terms of momenta 
Pi t 2 and metric tensors 

•■•Pi» k ---P2 H ---9w h --- ■ (39) 

Then, as (^3-^1,2) = 0, all tensor structures containing pi^ k or p 2flj will vanish when contracted 
with £% k or . Analogously, <7 WWl cancels when contracted with £% e £% h because l\ = 0. What 
proves the theorem. In the same way it can be shown the identity for t±. 



The third index is in parentheses to remind that the equation is valid for rank 2 and 3. 
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Corollary: 



/ ^zvW 9 • * = / ^W^" ' <W 2 * = ' (40) 

It can be proved again performing an explicit tensor decomposition. This theorem allows to 
replace in the integrand of a rank 2 or 3 three-point integral 

Q,Qu (q P ) - P [D X q a ] {9xat, u } (q p ) - 7 q 2 {t pu } (q p ) . (41) 

Indeed, Q P Q U in Eq. (J2TJ) read 

QpQv = (q ■ + (l ■ + (q ■ h) (q • 4) V ■ ( 42 ) 

Then, the previous theorem guarantees that the first two terms give zero after integration, 
and the corollary that the same remains true when they are multiplied by one and only one 
additional integration momentum q p . Finally, Eq. (|41|) and the same steps as in Sections 121 
and 01 result in the recursion relation ([35)1 . 

To conclude, we observe that similar methods can be also used to compute two-point 
tensors, as explicitly shown, for a particular case, at the beginning of next Section. However, 
the case m = 1 is so simple that we do not find any advantage with respect to standard 
reduction techniques. 



5 Rank one tensors 

As already observed at the end of Section Eq. cannot be applied, as it stands, to 
reduce tensors of rank one. In this Section we show how to cope with this situation. 



5.1 The m = 1 case 

The standard Passarino-Veltman decomposition gives 

4ri = f- 2 {ho i{ n) + iir\i) - ^ n) (o)} . (43) 

The same result can be derived extending our method. We can write 

Pi = *i + j with i\ 2 = and (h -£ 2 ) = 2 (£, ■ Pl ) = (£ 2 ■ Pl ) = p\ , (44) 

what corresponds to the second case of Footnote El with oti = 1/2. The massless 4-vectors 
£ h2 in Eq. (gU) always exist. For example, if p\ > 0, l x = (M/2, -M/2, 0, 0) and £ 2 = 
(M, M, 0, 0) in the frame where pi = (M, 0). The corresponding 4-vectors £3^ are defined as 
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in Eq. (fTHjh As besides 1^ must be proportional to Eqs. (fTHjl and (jjlj) imply 



JJqJJi 

Using them and Eq. (fTHjl one obtains 



d n q WW (q-£2) 



1 /> ? 2(g. Pl )rf 



D D 1 



and Eq. (g3J). 



5.2 The m = 2 case 

Using Eq. (fTHj) and the Theorem in Eq. (|37|) one can show that 

q* 1 P r „ -D M 



7 ■/ A)A 



so that 



r (n) _ x(n) 



where J is defined in Eq. f)l(J|l . 

5.3 The m = 3 case 

We use again Eq. (fTHj) to write 



/■(«) 



7 



27 



D Q D x D 2 Dz 

Multiplying and dividing (q ■ ^ 3)4 ) by (p 3 ■ £ 4i3 ) we can express 

Tr^^^l 



(?'4,4) 



(P3 ■ 4,3) 



(45) 



(46) 



(47) 



(4£ 



(49) 



(50) 



where the trace containing 75 vanishes upon integration because it is proportional to the 
totally antisymmetric tensor e(£x, q, £2^3) ■ Then, the following substitution is allowed in 
the integrand of Eq. (}4T?|) (where we also use Eq. (|TS|l ) 



-3A 



(P3 • ^4,3) 



(51) 
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from which the desired result follows (using Eqs. (fTUj) and (HHJ)) 



An) 
J 3;At 



P An) , 
7 



j30-<3 + -'2 



(P3 

(3) 



+ 



4 n) (o) 



2 ^„A An) 
7 ^ 3; A 



(52) 



5.4 The m > 3 case 

Contracting Eq. (JT8|l with p 3j4 one can write function of (q ■ Pi), i — 1, 



4: 



% = + ^ 



^ (pg (2 (g ■ p 4 ) - ^ m Z^ - ^4 f 2 (g • ft) - ^ft^ 



7 



7 



5 = (h-p4)(h-P3)-(h-P3)(h-P4). 

Then, after reconstructing the denominators and integrating, one gets 



(53) 



An) 
m; u 



7 ' h 2d 



- p2 



/ 3 o4 n) + &(3) - eii(o) - ^psx 4: ] 



\ P3 



(«) 



/ 4 o4 n) + 4i(4) - 7^(0) - ^ m JW 



2f3 



7 



2/9 



7 



(54) 



Note that \S\ 2 is proportional to the Gram determinant of the 4-momenta £i,i 2 ,P3,P4- The 
appearance of inverse square roots of Gram determinants is a peculiarity of our formalism 
and it will be discussed at length in Sectional 

We close this Section by observing that Eq. (J53J) could be used instead of Eq. (|2~4"j) when 
m > 3 to derive a recursion relation alternative to Eq. @. In fact nothing prevents from 
multiplying by an arbitrary number of 4-vectors q p - • -q T before integrating. The reason why 
we prefer Eq. (J24J) is because it involves only three out of the m external 4-momenta, while 
a fourth momentum is necessary to write down Eq. (|53|). This has important consequences 
when studying collinear or coplanar configurations, as we shall see in next Section. 



6 Study of exceptional configurations 

We are interested in the behaviour of our formulae at the edges of the phase-space, where 
two or more momenta can become nearly linearly dependent. We shall first show that only 
square roots of the two Gram determinants 

A = -A12 = (pi ■ p 2 ) 2 - p\p\ 
A 123 = 2(£ 1 -£ 2 )(£ 1 -p 3 )(£ 2 -p 3 )-p 2 3 (£ 1 -£ 2 ) 2 (55) 

appear in the denominators of Eq. ©. Furthermore, the only occurrence of a square root 
of a rank four Gram determinant is in Eq. (}53|) . These two facts make our approach numer- 
ically more stable than conventional methods, in which the scalar coefficients of the tensor 



12 



decomposition develop poles proportional to I/A12, I/A123 or 1/|<5| 2 at each step of the 
reduction. Secondly, we shall argue that Eq. (^J) is such that the numerical cancellations 
occurring among the tensor loop functions in the limit of exceptional momenta are kept as 
local as possible. Otherwise stated, each of the three terms of Eq. (jHJ) is separately well 
behaved when A12 — > or A123 — > 0. Although this cannot solve by itself all problems of 
numerical inaccuracy, it helps in decreasing the values of A 12 and A 12 3 for which approxima- 
tions to Eq. Q should be used. Thirdly, we shall show how to deal with configurations with 
exactly zero Gram determinants and give general prescriptions on how to cure the numerical 
instabilities occurring near the zeros of A12 and A123. 

We should then investigate all possible denominators appearing in Eq. Q in the limit 
of exceptional momenta. We start considering the 4- vectors in Eq. ()14|) . Inserting a 12 
in the definition of (3, one finds 

" = ± A< (56) 

where the sign depends on the sign in a\ (see Footnote^. Despite of this, remain well 
behaved in the limit A12 — > 0. In fact, splitting the external momentum p 2 as follows 

P2 = VPi + (P2 - VPi ) = VPi + <M > (57) 
and choosing rj and such that (p± ■ h) = and h 2 = —pl, one gets 



vA , , 

P2 = T)Pl H o - 71 ' ( 58 ) 

Pl 



with 



(Pl ' P2 



and n M = — = 
Pi VA 



P1P2 - (Pi -P2)P\ ■ (59) 



Therefore, in terms of p\ and n, we get well defined expressions even when A 12 — > 0: 

h = l(piTn), £ 2 = ^-( P i±n). (60) 

If needed, one can use higher numerical accuracy just in the computation of n M , or choose 
a particular frame to stabilise the result. For example, for time-like pi, one takes = 
M (0,p 2 /\p2\) when p^ = (M, 0). Finally, when both p\ vanish, a ii2 = as noticed in 
Footnote and l\p = pi >2 . When only one p\ is zero, £ 1>2 are still well defined provided 
(Pi " P2) 7^ 0, but this can also be cured, as discussed later. 

Let us now analyse the terms (3, I/7, 1/(^3 -^3) and 1/(^4 -^3) in Eq. Q in turn, where 
the last two quantities are hidden in the definition of the rank three tensor T^x- We start 
considering the case pf 7^ 0. Then (3 is proportional to 1/y/A, but the combination of tensor 
loop functions multiplying (3 in Eq. (jHJ) is always such that the product is well behaved in 
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the limit A12 — > 0. The proof is simple and follows from the fact that these terms come from 
defined in Eq. (fTHj). Indeed, it can be shown using Eq. (|50|l that 



D» = T 2^- [(q ■ n)W - (q ■ Pl )rf] , (61) 
Pi 



thus compensating the 1/vA pole coming from (3. Of course, we cannot simplify this pole 
analytically, because by doing so we cannot express D^ L back in terms of denominators. 
However, Eq. (|61J1 shows that well behaved combinations of loop functions naturally arise 
in our method. This is in contrast with the well behaved groupings of Ref. ^7], result of 
a very complicated compensation of logarithms and di-logarithms in groupings obtained by 
differentiating with respect to external parameters or by developing scalar integrals in 6 + e 
or higher dimensions. Obviously, such a compensation has also to occur after integrating 
over d A q. The nice feature of our approach is that this works before performing the actual 
integration, what keeps things much simpler. In addition, in our case the cancellations in 
the numerator should compensate a factor 1/y/A, and not 1/A. Finally, when at least one 
Pi is massless, j3 is simply 1. 

Next we concentrate on the zeros of 7. Since 

7 = Vlvl (62) 

(Pi • P2) =F V A 

there is no problem for massive P \ y2 . On the other hand, when pi >2 are both massless, 7 = 
2(pi ■ P2) can vanish. However, such configurations correspond to true collinear singularities 
of the amplitude, which are cut away in physical observables. When p\ — and p\ 7^ (jp\ 7^ 
and p\ = 0), 7 can become zero at the edges of the phase-space because 7 = 2(pi ■ p 2 ) = 
is a zero of the Gram determinant A pi _ P2iP2 (A pij p 2 _ pi ). For loop functions with m > 3, one 
simply renames P i ; s assigning the massless 4- momentum to p^. On the other hand, if m = 2 
and to fix things pi = 0, one redefines £^2 in Eq. (fTHjl using pi — p 2 and p 2 instead of pi 
and p 2 in order to move this pole to (3. The latter is important because as we have seen, f3 
behaves as the inverse square root of a Gram determinant, while poles in I/7 2 are present 
everywhere in our formulae. The net effect of this choice for £ 12 is the replacement 

^_ 1;p ... T (l) - ^- 1;p ... T (l) - /™-i ;P ... T (2) (63) 

in Eq. ©. 

We turn to 1/(^3 ■ pz) and 1/(^4 • P3). They are both proportional to l/-y/A 12 3. In fact, 
1(^4 -p 3 ) | 2 = |(4 -p 3 ) I 2 = v{t x )feuj- u{e 2 )v{i 2 )l)zoo- u(lx) = v{i x )falifau)-u{h) 

= ^{^hhU = (£-^) A i23- (64) 

This completes the proof of that only square roots of Gram determinants appear at each 
step of our reduction procedure. 
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Now, we shall show that also the last two terms of Eq. (^J) contain only well behaved 
combinations of loop functions in the limits A 12 — > and A 123 — ■> 0. To prove this, we split 
the 4- vector p 3 in a similar way as we did with p 2 in Eq. (|57j) 



P3 



Choosing rji, r\ 2 and 



mr 



such that {£\ ■ rh 



(£2 • J7i) — and m 2 = — (i 



(65) 



1 • £ 2 ) , one gets 



(£1 ■ l. 



V2 



(£l ■ £2 



i(£i-£2)P3 fl -(£2-p 3 )£S-(£i-P3)£2»] 



123 



Then 



(66) 



(67) 



(4,4 • Pa) = (4,4 • 77l) , 

implying that T^ u , defined in Eq. (J22J) and depending only on the ratio between (£ 3 -p 3 ) and 
(£4 'P3)) behaves smoothly in the limit A 12 3 — > 0. This tells us that the only singularity in 
the second term of Eq. (jUJ) can come from I/7, but this can be cured as explained above. 
To conclude, the last term of Eq. Q is proportional to T^ u \ ~ ^a 123 ' Expressing £i j2 in Eq. 
f!65|) in terms of pi t2 one obtains 



P3 





Pi + 


'2(3 ' 


— [r 2 -Pz) 


— {n -P3) 


. 7 







Then, using Eq. (JB'Hj) the coefficient of q x T^ u \ in Eq. 



V A123 * 

can be written as 



(6£ 



1 / - x V Al23 



(69) 



Therefore, the last combination of loop functions in Eq. Q must combine in such a way 
that the pole 1/a/A 12 3 coming from the tensor gets compensated. 

When A 12 or A 12 3 are exactly zero one has to rely on a different strategy. Let us con- 
centrate on the case when, for example, p 2 = Xpi exactly. If m is large enough, one simply 
chooses within the set {pj} a different subset of three independent 4-momenta to perform 
the reduction. D 2 acts then as a spectator denominator. When due to cancellations of 
denominators one is left with tensors integrals 



d n a - g/ i"- g "_ 
D D l D 2 D 3 ' 

one switches to the reduction valid for 3-point functions. In fact, when p 2 = Xp 1 , 



I 



d n q- q »'" qv 



D Q D 1 D 2 D Z 



and 



d n q q J l '-' q - 
q D D l D 3 



(70) 



(71) 
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share the same tensor basis built up in terms of and metric tensors. At the end, only 
2-point like tensors remain to be further reduced with the standard techniques of Ref. [2]. 
The same procedure also works when three or more momenta become linearly dependent. 
With the outlined method, reliable approximations can be easily obtained near the zeros 
of the Gram determinants, where numerical cancellations occur among the loop functions 
of each well behaved combination. When, for example, A 12 <C 1, one simply starts the 
reduction using D 2 = (q + Api) 2 — m 2 2 . This is completely equivalent to the Taylor expansion 
performed in Ref. [^J to extract the constant terms in each grouping. The advantage of 
our approach is that we do not need to explicitly develop tensor functions. Just at the end, 
when everything is reduced to scalar (4 + e) -dimensional loop integrals, only these need to be 
computed precisely, also in the regime of vanishing Gram determinants. This can be done, 
for example, as explained in Refs. |U |HJ EJ UZj • 

A problem which remains to be solved is the determination of the maximal values of 
A 12 , A 12 3 and S in Eq. (|55|l below which the approximations to Eq. (JSJ) should be used. As 
in Ref. [T7j such values can be only found performing dedicated numerical studies, and it 
is then difficult to give general prescriptions. We do not want to get deeply involved into 
the subject here, but just mention that by taking advantage of the fact that our reduction 
procedure takes place before integration, tests on the numerical stability of the formalism 
are possible without even evaluating the loop integrals. For example, for any given arbitrary 
4- vector q the integrands on the r.h.s. of Eq. (jHJ) should add up in such a way that at the 
end of the recursive algorithm the result is numerically equivalent to 

DqDi ■ ■ ■ D m 1 ' 

We have performed such a check on tensor integrals up to rank 4. 



7 Summary 

We have presented a method to compute numerically and recursively tensor integrals appear- 
ing in one-loop calculations, and relevant for the next generation of pp and e + e~ colliders. 
The treatment is applicable irrespective of the number of external legs to any configuration 
of internal and/or external variables, and only requires the knowledge of the standard set of 
scalar one- loop integrals. We distinguish the cases of 3-point tensor integrals (Section HJ), as 
well as of rank 1 (Section |3J), which are treated separately with similar techniques to those 
used in the general case (Sections |2] and |HJ). Singular kinematical configurations are analised 
in detail (Sectional), finding a smoother behaviour than in other approaches. In addition, 
we have studied all possible sources of numerical instabilities, giving general prescriptions 
on how to cure them. A code implementing the proposed method will be made available in 
the near future. 
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A The vectors ^ 3 4 

The 4-vectors £3,4 defined in Eq. (1TKJ) enjoy useful properties. By using the Dirac equation 
and the completeness relations for massless spinors one immediately derives 

(4,4 A2) = 0, 

1 



1 

-2 



■ v(£ x ) ^ u~ u{£ 2 ) v(£ 2 ) 7^ u~ u(£ 2 ) 

■v{£ x )^^-u{£ 2 ) = 0, 



v(£ 2 )Mh) 

£\ = v(£ 2 )^u-u(£ x )v(£ 2 ) lfl u>-u(£ x ) = 0, (73) 
where b is an arbitrary 4- vector different from £ Xt2 . Furthermore, 

(4- h) = v{£ 1 ) llx u-u{i 2 )v{l 2 )^u-u{t 1 ) = v{£ 1 ) 1 J 2l ^-u{£ l ) 

= ^Tr[f l7 ^27 M ] = -4(^-4), 

= ^[hiM = 4(q-£i)(q-£ 2 )-2q 2 (£ 1 -£ 2 ) } 
(q-h)(q-£ 3 ) = j^MJi^Uiiu-uih) 

2 {iq 2 (h-£2)-2(q-£ 1 )(q-£ 2 )}(b-£ 3 ) 



(b ■ Q 

+2 [(q ■ £ x )(b ■ £ 2 ) - (q ■ b){t x ■ £ 2 ) + (q ■ £ 2 ){£ x ■ b)](q ■ £,)} . (74) 

In the same way one computes (q ■ £±){q ■ £4). 

For numerical applications one needs £ 3i4 in terms of £ x>2 . Given £f = (£ i0 , £ ix , £ iy , £ iz ) and 
using |23j : 
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7 



v{£i) u 









' 


-a k \ 




o") 






o 1 



(fe^Cj ,0,0) , u(^) 



the desired result follows 





V 4 j 



7 



(7 





-a 



+ 



30 


= 6i6 2 


+ c{4 , 


4o 


= b 2 h 


+ c 2 cj , 


3a; 


= hc 2 


+ c[b 2 , 




= b 2 ct 


+ c 2 bi , 


32/ 


= «(cr 




£^y 


= i{ c 2 


b\ - b 2 cf) , 


3z 


= b x b 2 


~ c l c 2 1 




= b 2 bi 


— C 2 Ci . 



ni 



(75) 



(76) 



B Extra integrals 

In this Appendix we compute the extra integrals appearing in the n-dimensional version of 
the proposed reduction method 



r(n;2e) 



d n qq 



~2£ 1m ' ' ' QfJ-23 



Dn---D r 



(77) 



where £ > and 2s are non-negative integers. It is convenient to define a new index 
d=£+s+l— m to classify the integrals according to their dimensionality in powers of 
[mass] 2 . One can convince him/herself that tensor integrals have a non-zero contribution at 
(9(1), which is the order we are interested in, only when d > 0. For example, 



/■(»»; 2) 

2;^ 



ITT 

T 



(p^+P2^ + 0(e). 



(7S 



Their general expression can be easily written for d < 0: 

.2 ffMr--M 2s r(m-s-l) 



r(n; 2(m-l+d-s)) 



-27T 







r^ + OW whe nc / = 0, 
when d < , 



(79) 

but for 4-dimensional metric 



where the symmetric combination g M1 ... M2s is as in Eq. 
tensors. 

Let us discuss the scalar case (s — 0) in more detail. Decomposing the integration |2Sj 



J d n q = J d*qd € (i (<f = -fi 2 



(80) 



and after using Feynman parametrization, one gets 



j(n; 2(m-l+d)) 



) 2m+d m 2 



7T 



e/2 



r (f /2) r l" l - 1 + <i+ 2 lr 



-d- 



x / [da] m x!n + " 2) , 



with 



/POO m 
[da] m = / da ■ ■ ■ da m 8(1 - V) oc k ) , Af m = + M, 

fc=0 



m 



53 akPk , = a mo + ^ «fc( m fc - P*) • (82) 

k=l k=l 



We explicitly compute here one- loop integrals coming from tensors with rank at most 
equal to the number of denominators. This is enough for most practical calculations and 
gives rise to only one possibility with d — 1, namely m — 1. A straightforward calculation 
for this integral, gives 



2 , 2 Pi 

m i + m o~Y 



+ 0{e) . (83) 



When m,Q = ml — p\ — the integral vanishes in dimensional regularization. 

Furthermore, all scalar integrals with d = are easily evaluated by observing that 
J[da] m = l/T(m + 1): 

e ;2(m - 1)} = -™ 2 ^r^ + °( e ) when rf = • ( 84 ) 
I [m + 1) 

When d = -1, Eq. JSTJ) gives 

j0» ; 2(m-2)) = _ i7r 2^_p/ _ 2 + £A r r i _£ 

r(e/2) V 2/ V 2 

x | [da] m ^ 1+f) . (85) 

Since the constraint £ > in Eq. (J77J) implies m — 2 > 0, one could conclude that 

4 n;2(m " 2)) =0 + C(e) (86) 

i/ and on/?/ if the integral over the Feynman parameters in the second line of (185)1 does not 
contains poles in 1/e. One may wonder whether this can actually occur for there may be 
infrared and collinear divergences. A simple reasoning shows that this is never the case. 
In fact, when I > 0, the presence of q 2e in the numerator always raises the powers of the 
quadratic form X m in the Feynman parameter integrals with respect to the "standard" i = 
loop functions, therefore forbidding the presence of collinear or soft divergences. The same 
reasoning shows that 

j(n;2( m -l + d)) = Q + ^ when d<Q ( g7 ) 

However, it is extremely instructive to compute Eq. (|85j) directly. For the most collinear 
and infrared divergent fully massless box diagram in Fig. iBlone gets the finite expression 



/ 



[da} 3 — = — Ms/t) Hs/t) + Li 2 (l + t/s) + Li 2 (l + s/t)} . 

S + t 
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Pi 







(P3 -Pi) 2 = (P3-P2)' 

s=Pl, f = (Pi - P2) 2 



P3 "Pi 



P3 -P2 

Fig. [Bj Fully massless box diagram. 



Since all higher point loop functions can be rewritten as combinations of box diagrams jlj [7] , 
this result explicitly proves Eq. (|57ji. 

We close this Appendix noting that, in any case, a contribution 0(1) can only develop 
for non- negative powers of the quadratic form X m , so that the integrand in the Feynman 
parameter integral is always polynomial. 



C Three-point tensors: general case 

This Appendix extends Eq. (|35p to three-point tensor integrals of rank higher than 3. We 
redefine 

dp = and Lii = --^Qn, (89) 

so that the first line of Eq. (jl8j) simply reads = + L M . Then 

i i 

9a»i • • • %i = E • • ■ L Vk-i d ^ q^ k+1 •••?«}+ II L f * • ( 90 ) 

k=l fc=l 

The first term always contains reconstructed denominators and we do not elaborate it any 
further. Inserting the definition of in Eq. (JTSj) in the second term this can be written 

fl W k = (-^Ytiv ■ h)\q ■ Uf~ k) Sl.. m , (91) 
fc=l \ z '/ fc=0 

where the tensor S*,.^ is defined to be the sum of all possible tensor products of i 4- vectors 
£3 and £4, such that £4, appears k times. For example, 

*-Vi/^2At3 = ^4/ii^i 2 ^3At 3 + ^4/ii^3At2^4/i 3 + ^4/i 2 ^4/i 3 • (92) 

The two terms with fc = and = i in Eq. (}9~T| do not contribute to the integral, due to 
the Theorem in Eq. (p?7|) . The product (g • £s)(g ■ £4) = /3q x D\ — jq 2 can be factorised out 
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from all the remaining terms, giving as result 



s 



k=l 



ot\---a.i-2\ Mi'-'Mi 



— I [{3q*D x - ig 2 } <T ■ ■ ■ q a ^ S [ai ... ai _ 2] 



t-i 



Eoi— k— 1 oO 



(93) 



fe=i 



Inserting Eq. (}9*3*j) into Eq. ([90)1 . dividing by the three denominators and integrating over 
d n q we obtain the desired relation 



An) 

'2; //v/Lti 



X 



7 fc=l 
(o j(n)\ 



27, 



•«i-2 



-7 



2 r(n) ,(n) /pv _ j{n;2) 11 



(94) 



where S^.'.V^* -2 ] is given in Eq. ()93|). but with the a indices lowered, and as in Eq. (J33|) the 
tensors in the numerator of the n-dimensional integrals are 4-dimensional. 
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